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Dynamics of crater formations in immersed granular materials 
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We report the formation of a crater at the free surface of an immersed granular bed, locally 
crossed by an ascending gas flow. In two dimensions, the crater consists of two piles which develop 
around the location of the gas emission. We observe that the typical size of the crater increases 
logarithmically with time, independently of the gas emission dynamics. We describe the related 
granular flows and give an account of the influence of the experimental parameters, especially of the 
grain size and of the gas flow. 



PACS numbers: 83.80.Fg : Granular materials, rheology, 47.57.Gc 
: Hydrodynamics, applied fluid mechanics 



Granular flow, complex fluid, 47.85.Dh 



I. INTRODUCTION 

Craters are part of the widespread phenomena ob- 
served in nature, going from impact meteorite craters 
to volcanic structures. Studies of crater morphologies 
have a wide range of applications, going from puzzling 
crater formation in drying paint [l[ to molecular dynam- 
ics @, 0, 0|- Among the main applications to natural 
phenomena, aside from meteorite impact crater, are the 
formation and growth of volcanic edifices, by successive 
ejecta emplacement and/or erosion. The time evolution 
and dynamics play a crucial role here, as the competition 
between volcanic-jet mass-flux (degassing and ejecta) and 
crater-size evolution may control directly the eruptive 
regime 5] . Attempts have been made to model the talus 
development of volcanic caldera through erosion [6j, or 
to constrain the morphology and dynamics of pyroclastic 
constructs via granular-heap drainage laboratory exper- 
iments 7]. 

Crater morphology in dry granular material has been 
extensively studied, both experimentally and theoreti- 
cally d 1, M, El E G3- Most of these studies in- 
vestigate the final, steady, crater shape resulting from 
the collision of solid bodies with the material surface and 
scaling laws are derived [Oj]. Note however that some 
authors also reported experimental study of the crater 
formation dynamics, including growth and collapse af- 
ter impact, in the "single impacting body" configuration 
[TEl ITq ] . In a recent work, Wu et al. have extended these 
studies to a particles stream impacting a dry granular 
bed [13 • 

In immersed granular material, one reports craters 
generated by an underwater vortex ring [18j . involving 
fluidized ejecta dynamics, or underwater impact craters 
generated by landslide [19(. Craters in immersed granu- 
lar materials can result either from two-phase or three- 
phase flows. In particular, water or gas flowing through 
an immersed granular bed can induce localized instabil- 
ities and fluidization [2(3, HH , which eventually leads to 
the formation of craters at the free surface. 

In a previous experimental study, Gostiaux et al. [22| 
have investigated the dynamics of air flowing through an 
immersed granular layer. They reported that, depending 



on the flow rate, the system exhibits two qualitatively dif- 
ferent regimes: At small flow rate, the bubbling regime 
during which bubbles escape the granular layer quite in- 
dependently one from another; At large flow rate, the 
open-channel regime which corresponds to the formation 
of a channel crossing the whole thickness of the granu- 
lar bed through which air escapes almost continuously. 
At intermediate flow rate, a spontaneous alternation be- 
tween these two regimes is observed. Interestingly, they 
noticed the appearance of a crater around the locus of 
air release but did not provide any extensive study of 
its formation dynamics. Here, we extend these seminal 
observations to a quantitative study of the resulting de- 
formation of the free surface of the granular bed. To do 
so, we reproduce the previous experimental conditions in 
two dimensions: In a vertical Hele-Shaw cell, the crater 
is then formed by two sand piles which grow and move 
away from each other as time passes by. By monitoring 
the evolution ot the free surface through time, we inves- 
tigate the effect of the different gas flow regimes on the 
crater dynamics. The results point out a grain-transport 
mechanism that differs significantly from the one involved 
in dry-sand dune motion [23|, [24], Hfl [26| . 



II. EXPERIMENTAL SETUP AND PROTOCOL 

The principle of the experiment is to force air, injected 
locally at the base of an immersed granular bed, to flow 
vertically across the material. The experimental setup 
thus consists in a 2-mm-thin vertical-cell containing the 
grains and water (Fig. [J) . The two vertical walls (glass 
plates) are 40 cm wide and 30 cm high. 

The air injection is assured by a system similar to that 
used in [22j: Air is injected, at a constant flow rate, in a 
chamber connected to the bottom of the cell. In practice, 
a small hole (diameter 1 mm) connects, at the center, the 
bottom of the vertical cell to the chamber, partially filled 
with water in order to tune the inner gaseous volume 
V. The air flow is controlled thanks to a reducing valve 
which maintains a constant overpressure at the entrance 
of a capillary tube: Provided that the pressure drop in 
the capillary tube is much larger than the fluctuations 
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FIG. 1: Experimental setup: Air is injected at constant flow rate $ into a chamber (volume V) connected to a thin vertical 
cell containing an immersed, initially flat, granular layer. Through time, symmetrical piles form, grow and move in opposite 
direction on either side of the crater. The experiment consists in observing the deformation of the granular-bed free-surface 
and in monitoring the overpressure AP inside the chamber, in order to account for the influence of the degassing process on 
the crater dynamics [the picture here illustrates the bubbling regime (see text).] 



of the pressure inside the chamber, the air flow toward 
the chamber is not significantly altered by the dynamics 
of the granular layer and remains constant. The flow 
rate $ (ranging from 1.5 mL/s to 4.3 mL/s) is tuned by 
changing either the capillary tube (rough) or the pressure 
difference imposed by the reducing valve (fine) and is 
subsequently measured (to within 0.1 mL/s) by means of 
a home-made flow meter. A sensor (MKS Instruments, 
223 BD-00010 AB), connected to an acquisition board 
(National Instruments, PCI-6251), monitors the pressure 
difference, AP (to within 10 Pa), between the chamber 
and the outside atmosphere. 

In order to account for the dynamics of the free surface, 
we image the system from the side. In a first configura- 
tion, a transparency flat viewer (Just NormLicht, Classic 
Line) located behind the cell is used to achieve a homoge- 
neous lighting of the whole bed which is imaged with the 
help of a Webcam (Logitech, QuickCam Express) con- 
nected to a PC. A small software (Astra Image Webcam 
Video Grabber) makes it possible to take one image of 
the system every 10 seconds and, thus, to record the dy- 
namics during several hours (typically 24 hours). In a 
second configuration, in order to observe the motion of 
the grains, we light up the sample from top with a linear 
light source (Polytec, DCR3) and take pictures with a 
high- resolution digital-camera (Nikon, D200). 

Initially, <E> is set to zero. Grains and water are intro- 
duced in the cell. The grains consist of spherical glass- 
beads (USF Matrasur, sodosilicate glass) that are sieved 
in order to control their size (typical diameters, 100-125 
or 400-500 /xm). Initially, after sedimentation, the grains 
sit at the bottom and we make use of a thin rod to level 
the immersed granular layer. We denote h g the height of 
grains above the bottom of the container and h w , that of 



the water free surface above the granular bed. 

The initial condition consists thus in a horizontal and 
flat layer of grains. After the opening of the valve, one ob- 
serves that air is creating several paths between the hole 
at the bottom and the free surface of the granular bed. 
During this transient regime, some air channels merge, 
or cannot reach the upper layer of grains (In this case, 
air bubbles remain trapped within the granular layer.) 
After several minutes, air crosses the granular bed along 
the vertical and bubbles are emitted in the water from 
an almost fixed position. Then, through time, a crater 
forms, the two piles on both sides growing and moving 
apart one from the other (see Fig. [1} . 



III. RESULTS 



We shall report the dynamics of the crater formation 
in regard to the gas flow regimes. First, we shall describe 
qualitatively the gas flow regimes (Sec. IIII A 1[) and the 
associated mechanisms of the gas emission at the free 
surface (Sec. IIII A 21) . Then, we shall describe qualita- 
tively the crater formation (Sec. IIII A 3[) and make use 
of these first observations to estimate the profile of the 
grain deposition. We shall also discuss the effects of the 
finite depth of the granular bed and of the finite water 
height (Sec. IIII A 4|) . We finally quantitatively study the 
influence of the gas flux (Sec. IIII B ip and of the grain 
size (Sec. IIII B 21) . 
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A. Qualitative observations. 

1. Gas flow regimes 

As already reported [22j], depending on the flow rate <E>, 
two main regimes of the air flowing through the granu- 
lar layer are observed. (The existence of the two regimes 
and the transition between them are tightly related to 
the non-Newtonian rheology of the granular material as 
proven by an experimental study conducted in another 
non-Newtonian material 27}.) On the one hand, the bub- 
bling regime, which is typically observed at small <f>, is 
characterized by a regular emission of successive bub- 
bles, independent from one another. In this regime, the 
pressure signal exhibits successive rises and drops, the 
latter being associated with bubble emissions from the 
injector at the base of the granular layer. On the other 
hand, in the open-channel regime, which is typically ob- 
served at large 4>, the system sustains a continuous air- 
flow through a channel crossing the whole granular layer. 
The overpressure AP associated with this continuous air 
emission is almost constant. In an intermediate range of 
one observes a spontaneous alternation between the 
two regimes: The channel forms after the emission of sev- 
eral bubbles and subsequently spontaneously pinches off 
after a finite time, leading the system back to the bub- 
bling regime. As a consequence, activity and rest periods 
are observed in the pressure signal [Fig. ISJa)]. The phe- 
nomenon is explained by the ability of the material to 
sustain a stable channel thanks to its peculiar rheology. 
At last, we point out that the deformation of the free 
surface of the granular bed does not seem to alter the 
gas flow process: It has been previously shown that the 
gas emission is mainly governed by local events occuring 
close to the free surface 22] . 




FIG. 2: Overview on the gas flow process: (a) At in- 
termediate flow-rate $, the pressure signal exhibits a spon- 
taneous alternation between activity and rest periods which 
correspond respectively to the bubbling and open-channel 
regimes. The emission of the bubbles at the free surface ac- 
cordingly exhibits two qualitatively different behaviors : (b) 
In the bubbling regime, a large bubble forms underneath the 
free surface of the granular bed and then bursts. The explo- 
sion pushes out of the bed a large number of grains, which 
are subsequently advected upwards in the bubble wake, (c) In 
the open-channel regime, the walls of the channel are stable 
and a continuous gas flow escapes the granular bed. In this 
regime the grains are only torn out from the granular bed by 
the water flow and subsequently advected in the wake of the 



small bubbles that form in water (d 
h w = 16 cm, and $ = 2.5 mL/s). 



in large range of <f>. 



400 fim, h g 



3. Geometrical description of the crater 



2. Bubble emission mechanisms 

A close look at the free surface of the granular bed 
points out two different air-release mechanisms associ- 
ated with the bubbling and open-channel regimes, re- 
spectively. In the bubbling regime, a gas bubble, while 
growing underneath the free surface, pushes up a thin 
layer of grains which forms the bubble head. Once it 
has crossed the interface, the bubble, while it rises up in 
the water, advects the grains in its wake [Fig. Hfb)]. In 
the open-channel regime, the air is released continuously 
through the channel whose walls remain at rest. The 
grains are advected upwards, from the free surface, by 
the water flow behind the small ascending gas bubbles 
[Fig. [SJc)]. From these observations, one could wonder 
about the relative efficiency of the two regimes in forming 
the crater. In particular, one could expect the bubbling 
regime to lead to a faster growth because the explosive 
bursting of the bubble apparently lifts a larger quantity 
of grains. In Sec. IIII B[ in order to answer the question, 
we shall report a quantitative study of the crater growth 



The ejection of the grains from the free surface and 
the subsequent deposition of the granular ejecta lead to 
the rapid formation of a crater, which is formed by two 
granular piles symmetrically positioned on each flank. 
Due to the permanent grain transport, the crater grows: 
The piles height increases while they move symmetrically 
away from the center (Fig. [3]). 

Before reporting quantitative measurements of the 
crater dynamics, let us mention some qualitative features 
of the crater formation. First, let us denote a and (3 the 
maximum angles that the piles make with the horizon- 
tal, respectively, inside and outside the crater. From di- 
rect observation [Fig. [5] , one can notice immediately that 
a and (3 can differ significantly, especially at the early 
stages of the crater formation [Fig. |Ua)]. We note that 
the inner flanks of the crater are almost straight and that 
a (about 28°) remains constant during the whole crater 
growth. In contrast, the outer slopes are not straight, es- 
pecially at the early stages of the crater formation, and 
(3 evolves in time. 

This behavior can be qualitatively accounted for by 




FIG. 3: Temporal evolution of the crater: The black 
lines are the result of the automatic detection of the free- 
surface profile. The angles a and (3 denote, respectively, the 
maximum angle the inner and outer flanks make with the 
horizontal. We define L, the distance between the two peaks, 
and x, the distance from the center (d = 100 fim, h g = 8 cm, 
h w — 16 cm, and $ = 2.2 mL/s). 



considering the grain flows (Fig. [5]). Once lifted up by 
liquid flow, the grains are pushed away from the center 
by the liquid flow which mainly consists of two large con- 
vective rolls and then, subjected to gravity, deposit back 
onto the free surface at finite distance from the center. 
If the local slope is smaller than the angle of avalanche 
the grains do not move anymore once deposited. To the 
contrary, if the local slope is larger than the angle of 
avalanche, the grains flow downwards along the slope. 
At the center the grains are very locally torn off from 
the surface by the upward liquid flow. The local slope 
almost immediately exceeds the angle of avalanche and 
grains flow downwards on both inner flanks to replace 
the granular material missing at the center. Thus, if a 
grain deposits onto the inner flank, it flows downwards 
toward the center, which again explains why a remains 
almost constant and only fluctuates between the angle 
of repose and the angle of avalanche [Fig. BJa)]. I n con ~ 
strast, at the early stages of the crater formation, the 
outer flanks only result from the deposition of the gran- 
ular material. The local slope of the outer flank being 
everywhere smaller than the angle of avalanche, a grain, 
once deposited, does not move anymore. However, at a 
finite time t, (3 reaches the value of the angle of avalanche, 
which results in surface flows along the outer flank. As 
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FIG. 4: Temporal evolution of the crater geometry: (a) 

Inner and outer angles a and /3 vs time t (5> = 1.67 mL/s). 
The grey zone marks the early stages of the crater growth 
during which the crater size compares with the size of the 
bubbles. The angle a, which is then not measurable, remains 
constant during the whole experimental time. In contrast, the 
angle /3 increases continuously to reach a constant value after 
a finite time. In order to show that the slope is not constant 
along the outer flanks, we report the temporal evolution of 
the slope at x = \L{t) + 1.5 cm: the local slope remains 
smaller than /3 until reached by the avalanche flow. Lines 
are only guides for the eye. (b) Local height h(x, t) vs time 
t ($ = 2.93 mL/s). Far away from the center, the height h 
of the granular layer increases due to the deposition of the 
grains. We observe a linear increase of h(x, t) with time until 
the region is reached by the grains avalanching along the outer 
flank, (c) Local growth velocity dh/dt(x,t) vs distance x 
(<!> = 2.93 mL/s). The data are successfully interpolated by 
a-^- exp (— #-) for large x, we get at a rough estimate L c — 
(1.13 ± 0.04) cm. (d) Distance L between the summits vs 
time t ($ = 4.26 mL/s). Whereas the crater forms quite fast, 
we observe a drastic decrease of the growth velocity at large 
times. Inset : L vs \n(t). At large times, L increases almost 
logarithmically with time t (d — 400 /im, h g = 8 cm, and 
h w — 16 cm). 



a consequence, j3 evolves in time [Fig. [4ja)] : Starting 
from a small value at the early stages of the crater for- 
mation, [3 increases to reach a constant value after a finite 
time. We note here that f3 then remains slightly smaller 
than a, which is probably explained by the fact that a 
significantly exceeds the angle of avalanche due to the 
continuous grain flow along the inner flanks. 

At this point, it is particularly interesting to make use 
of these observations to estimate the spatial distribution 
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Colored grains 



FIG. 5: The granular flows: The second configuration is 
used: we light up the sample from top with the linear light 
source and take pictures with the high-resolution digital cam- 
era. Initially, a thin layer of colored grains is deposited at the 
free surface of the initially flat and horizontal bed. Then, a 
series of ascending rolls pushes the grains away from the verti- 
cal central axis. Subsequently, the grains gently deposit back 
onto the free surface of the bed. Along the piles flanks, pro- 
vided that the local angle exceeds the angle of avalanche, we 
observe continuous avalanches. Inside the crater, the flowing 
granular-material partly replaces the grains advected upwards 
at the center. Along the outer flanks, the deposited granu- 
lar material either flows or sits at the free surface. One can 
clearly observe, in the inset, that new material, deposited far 
away from the center forms a thin layer of material which re- 
mains at rest. The deposition flow rate, Q(x), is estimated 
from the temporal evolution of the growth velocity dh/dt of 
the deposited layer in that region [Fig. |4jc)] (d — 400 /im, 
h g = 8 cm, h w — 16 cm, and $ = 3.5 mL/s). 



of the deposited grains. To do so, let us now consider the 
local height h(x,t) of the free surface at the distance x 
from the center at time t [Fig.QJb)]. Close to the center, 
h(x, t) exhibits a complex behavior which results from 
both the deposition and the surface flows. Far from the 
center, h(x,t), which results from the deposition alone, 
evolves linearly in time. From this latter observation, we 
deduce that the advection is not significantly affected by 
the crater growth and that the distribution of the de- 
posited grains, far away from the center, is almost con- 
stant in time. At an intermediate distance from the cen- 
ter, we observe a change in the temporal evolution of 
h(x, t) : At small times, the dynamics is only due to the 
deposition whereas later, when f3 reaches the avalanche 
angle, the local dynamics results from both the deposi- 
tion and the surface flows. The volume, Q(x), of granular 
material deposited at the distance x from the center per 
unit distance and per unit time is proportional to the ini- 
tial slope dh(x, t)/dt far away from the center [Fig.fJJc)]. 
Seeking for a simple mathematical description of the de- 
position flow, we guess that the grains, pushed away from 
the center by the liquid flow, have a negligible probabil- 
ity to deposit back at the center and we propose to in- 



terpolate dh(x,t)/dt by a-^- exp (— -^-), which makes it 
possible to extract a characteristic length L c at which 
the grains are deposited away from the center. We shall 
show in Sec. IIII Bl that the shape of the distribution Q(x) 
accounts for the evolution of the distance L(t) between 
the two summits as a function of the time t [Fig. 0Jd)] . 



4- Finite size effects 

We aim at reporting the growth of the crater in an 
infinite surround, not limited by the finite depth of the 
granular bed, h g , or by the finite water height, h w . 

The growth of the crater stops when the thickness of 
the granular bed at the center vanishes. Taking into ac- 
count the angle of avalanche and the conservation of the 
granular-material volume, we can estimate the maximum 
accessible value of L, L 9 max = (2 + \/2)h g / tana. We 
shall report experimental results obtained for h„ = 8 cm 



so that L r . 



50 cm is about the lateral width of the 



cell. The growth of the crater shall thus not be limited 
by the finite depth of the granular bed. 
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FIG. 6: Distance L(30min) vs water height h w : For h w < 

L/4, the growth of the crater is limited by the water height 
and L(30min) increases linearly with h w . In contrast, for 
h w > L/4, the growth of the crater is not altered by the water 
height and L(30min) does not depend on h w (d — 100 /im, 
h g — 8 cm, and $ = 3.06 mL/s). 



In the same manner, the height of the piles can obvi- 
ously not exceed h w . We estimate, the maximum accessi- 
ble value of L, L^ ax = 2(1 + \f2)h w j tan a. However, we 
checked experimentally that L^ ax largely underestimates 
the finite water-height effect: Reporting L at a given, 
large, time t = 30 min as a function of h w , we obtain 
experimentally that the growth of the crater is limited 
by the water height for h w < L/4 (Fig. [6]). Interestingly, 
we observe, in addition, that 1,(30 min) does not depend 
on h w if h w > -L/4 : The crater growth is not altered by 
the water height if the latter is large enough. We shall 
report experimental results obtained for h w = 16 cm so 
that L^ ax ~ 64 cm is larger than the lateral width of the 
cell. The growth of the crater shall thus not be limited 
by the finite water height. 
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pend significantly on the grain size (the slope in the 
semilog plot remains almost unchanged) and we estimate 
Lq = (3.5 ± 1) cm for d = 400 /im and d = 100 /mi. In 
contrast, we estimate Vq ~ 0.4 mL for d — 400 (im and a 
significantly smaller value Vq ~ 0.02 mL for d = 100 /im. 
The precise dependence of Lq and Vq on d is difficult 
to access experimentally as such measurements require 
precise knowledge of the origin of time (a delay alters 
significantly the slope and the offset in the semilog plot) 
which we are missing. Nevertheless, Vq is observed to de- 
crease drastically with the grain size whereas Lq remains 
almost constant. 



FIG. 7: Distance L vs. time t: The distance L between the 
summits increases logarithmically with time, independently 
from the gas emission regimes. Indeed, at small <!>, the gas 
flow regime consists mainly in the emission of independent 
bubbles, whereas it consists mainly in a continuous emission 
of gas at large <&. From the interpolation of the experimental 
data with L/Lq = \n(&t/Vo), we estimate Lo — (3.5±0.7) cm 
and Vq = (0.4±0.2) mL for d = 400 /im (black curves) whereas 
we estimate Lo = (3.5 ± 1) cm and Vb = (0.02 ± 0.01) mL for 
d — 100 ^m (grey curve) (h g = 8 cm, h m = 16 cm) . 

B. Influence of the gas flow $ and grain size d 

1. Dependence on the air flux "!> 

We report the distance L between the piles as a func- 
tion of time t, for different air flow rates $ (Fig. [7]). We 
observe that L increases logarithmically with t and we 
write L/L = ln($t/Vb), Vb being a typical volume. 

In the bubbling regime, when the gas emission at the 
free surface consists of the periodic emission of indepen- 
dent bubbles at the free surface, a simple argument can 
account for the proposed dependance on the total emit- 
ted gas volume, &t, alone: The dimensions of the piles, 
length L or height, are expected to depend on the total 
number N g of grains displaced by the water flow. Thus, 
N g being proportional to the number of emitted bubbles, 
which itself, at a given time t, is proportional to &t, we 
expect L to be a function of <frt alone, as observed experi- 
mentally. In addition, we observe that the relation holds 
true even at large $ when the system exhibits, almost 
always, the open-channel regime. Thus, contrary to the 
intuition, the small bubbles emitted by the open-channel 
are as efficient as the exploding bubble of the bubbling 
regime in lifting the grains: This observation is again in 
agreement with the fact that the grains are mainly lifted 
by the water flow and not expelled by the explosion of 
the bubbles 

2. Dependence on the grain size d 

From the interpolation of the exprimental data, one 
can estimate the typical length Lq and the typical vol- 
ume Vq. One observes that the length Lq does not de- 



IV. DISCUSSION 

The crater growth is due to the advection of the grains 
and to the subsequent deposition away from the location 
of the gas emission. From the observation of the deposit 
far away from the center (Sec. IIII A3|) . we propose that 
the deposited flow can be written, for large distance x, 
Q(x) = a j- exp (—#-)• This relation accounts for the 
logarithmic growth of the crater. Indeed, assuming the 
angle a on both sides of the piles, a simple geometrical 
analysis gives the volume, v — \ j^^^L 2 of the grains 

on the outer side of one pile (x > L/2). Noting that only 
the deposition of grains on the outer flank of the pile con- 
tributes to its growth, we can write ^ — J^j 2 Q( x ) dx, 
which leads to the equation governing the pile growth, 

1 tana dL f°° x ( x \ 

— = — L— — = a — exp — — dx. (1) 

4(1 + V2) 2 dt J L/2 L c y \ L c ) K > 

Provided that L L c (which is almost always sat- 
isfied experimentally), we obtain that, asymptotically, 
L increases logarithmically with time t according to 
L{t) ~ L Q \n{$t/V ) with L = 2L C and $/Vq = 
a(l + v / 2) 2 /tana. The experimental rough estimates, 
obtai ned for d = 400 /mi, L c = (1.13 ± 0.04) cm 
(Sec. IIII A 3|) and Lq = (3.5 ± 0.7) cm (Sec. IIUBlj) are 
in fair agreement with this expectation. The same con- 
clusion holds true for the smaller grains, d — 100 /im, as 
we measured L c = (1.38 ± 0.04) cm for the same value of 

Finally, if we interpret the volume Vb as the gas volume 
necessary to lift or move a given quantity of grains, one 
can easily understand that Vb decreases with d. Indeed, 
the advection process, whose associated force scales like 
d, competes with the buoyancy force which scales like 
d 3 . As a result, large grains are more difficult to lift than 
small ones and the gaseous volume necessary to move 
them away is larger. 

The advection process remains difficult to model in de- 
tails. Indeed, the liquid flow in the wake of the bubbles 
is generally turbulent and it is barely possible to account 
for the advection of the solid particles in such a complex 
stream field. There also remains a very important open 
question: We observed that the growth of the crater is 
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not altered by the water depth (if large enough). How- 
ever, we observe that the vertical size of the large-scale 
convective rolls compares to h w , whereas their lateral size 
compares to L/2. One thus would expect the typical dis- 
tance L c to depend on h w and, even, on time. Our study 
does not display such an effect. The typical vertical dis- 
tance to be taken into account is probably the distance 
over which a grain follows the bubble in its wake, in other 
words, the maximum accessed altitude. 

In conclusion, we reported experimental data on the 
dynamics of a crater growth at the free surface of an 
immersed granular crossed by an upward gas flow. We 
observed that, due to the peculiar transport of the grains 
in the wake of the rising bubbles, the typical size of the 
crater increases logarithmically with time. The dynam- 



ics is demonstrated not to be altered by the gas flow 
regime and only to depend on the overall gas flow rate 
and the typical size of the grains. However, even if we 
clearly established that the global dynamics of the crater 
is compatible with the spatial structure of the grain depo- 
sition around the gas-emission locus, we are still missing 
a complete modeling of the grain advection by the tur- 
bulent flow, which deserves to be further investigated. 

ACKNOWLEDGMENT 

G.V. acknowledges a grant by CONICYT (Comision 
Nacional de Investigation Cientifica y Tecnologica, Gob- 
ierno de Chile). 



P. L. Evans, L. W. Schwartz, and R. V. Roy, J. Colloid 
Interface Sci. 227, 191 (2000). [15 
Z. Insepov, L. P. Allen, C. Santeufemio, K. S. Jones, and 
I. Yamada, Nucl. Instrum. Methods Phys. Res. B 206, [16 
846 (2003). 

Y. Yamaguchi, J. Gspann, and T. Inaba, Eur. Phys. J D [17 
24, 315 (2003). 

T. Aoki, and J. Matsuo, Nucl. Instrum. Methods Phys. [18 
Res. B 241, 594 (2005). 

A. W. Woods, and S. M. Bower, Earth Planet. Sci. Lett. [19 
131, 189 (1995). 

H. Obanawa, and Y. Matsukura, Geomorphology 97, 697 [20 
(2008). 

C. Riedel, G. G. J. Ernst, and M. Riley, J. Volcanol. [21 
Geotherm. Res. 127, 121 (2003). 

Y. Grasselli, and H. J. Herrmann, Gran. Matt. 3, 201 [22 
(2001). 

J. S. Uehara, M. A. Ambroso, R. P. Ojha, and D. J. [23 
Durian, Phys. Rev. Lett. 90, 194301 (2003). 

D. Lohse, R. Bergmann, R. Mikkelsen, C. Zeilstra, [24 
D. van der Meer, M. Versluis, K. van der Weele, M. 

van der Hoef, and H. Kuipers, Phys. Rev. Lett. 93, [25 
198003 (2004). 

X.-J. Zheng, Z.-T. Wang, and Z.-G. Qiu, Eur. Phys. J. E [26 
13, 321 (2004). 

K. Wada, H. Senshu, and T. Matsui, Icarus 180, 528 [27 

(2006) . 

S. J. de Vet, and J. R. de Bruyn, Phys. Rev. E 76, 041306 

(2007) . 

A. M. Walsh, K. E. Holloway, P. Habdas, and J. R. 



de Bruyn, Phys. Rev. Lett. 91, 104301 (2003). 

J.-F. Boudet, Y. Amarouchene, and H. Kellay, Phys. Rev. 

Lett. 96, 158001 (2006). 

S. Yamamoto, K. Wada, N. Okabe, and T. Matsui, Icarus 
183, 215 (2006). 

S. M. Wu, H. P Zhu, A. B. Yu, and P. Zulli, Gran. Matt. 
9, 7 (2007). 

A. Suzuki, I. Kumagai, Y. Nagata, K. Kurita, and O. S. 
Barnouin-Jha, Geophys. Res. Lett. 34, L05203 (2007). 
H. M. Fritz, W. H. Hager, and H.-E. Minor, Exp. Fluid 
35, 505 (2003). 

P. Rigord, A. Guarino, V. Vidal, and J.-C Geminard, 
Gran. Matt. DOI 10.1007/sl0035-005-0214-l (2005). 
F. Zoueshtiagh, and A. Merlen, Phys. Rev. E 75, 056313 
(2007). 

L. Gostiaux, H. Gayvallet, and J.-C Geminard, Gran. 
Matt. 4, 39 (2002). 

R. A. Bagnold. The Physics of Blown Sand and Desert 
Dunes. Chapman and Hall, London (1941). 
K. Kroy, G. Sauermann and H.J. Herrmann, Phys. Rev. 
Lett. 88, 054301 (2002). 

B. Andreotti, P. Claudin and S. Douady, Eur. Phys. J. B 
28, 321 (2002). 

P. Hersen, K.H. Andersen, H. Elbelrhiti, B. Andreotti, P. 
Claudin and S. Douady, Phys. Rev. E 69, 011304 (2004). 
T. Divoux, E. Bertin, V. Vidal and J.-C. Geminard, 
Intermittent outgassing through a non-Newtonian fluid, 
larXiv:0810 .3095vl [cond-mat.stat-mech] (2008). 



